SQUIRALS AND BEYOND: SUBSTITUTION TILINGS 
WITH SINGULAR CONTINUOUS SPECTRUM 



MICHAEL BAAKE AND UWE GRIMM 

Abstract. The squiral inflation rule is equivalent to a bijective block substitution rule and 
leads to an interesting lattice dynamical system under the action of Z 2 . In particular, its 
balanced version has purely singular continuous diffraction. The dynamical spectrum is of 
mixed type, with pure point and singular continuous components. We present a constructive 
proof that admits a generalisation to bijective block substitutions of trivial height on Z d . 



1. Introduction 

Dynamical systems of translation bounded measures on R rf with pure point spectrum are 
well understood by now. This owes a lot to the equivalence between pure point dynamical 
and diffraction spectra [331 E]) an d to the general characterisation of model sets via dy- 
namical systems; see [12J and references therein. Prom an applied perspective, pure point 
spectra are linked to the analysis of periodic and almost periodic systems, such as crystals 
and quasicrystals, by standard crystallographic methods; see [18] for background, and [TO] 
for a recent review in the aperiodic setting. In fact, the understanding of pure point systems 
has improved sufficiently that one can begin to attack the corresponding inverse problem 
systematically [35] . 

As soon as one enters the realm of mixed spectra, the picture is less transparent. While 
some partial understanding exists for systems with absolutely continuous spectral components 
(for instance through the connection with the highly developed field of stochastic processes; 
see [5] and references therein), one is pretty much in the dark when it comes to systems with 
singular continuous spectra. The classic Thue-Morse system, see [2] for background, and its 
generalisations in the spirit of [30\ [6] are notable exceptions, which are all one-dimensional. 

Beyond the theoretical interest, there was little motivation in the past to dive further into 
systems with continuous spectral components, at least not from an applied point of view. 
However, with the modern measurements possible in materials science, spectra of this type 
are detectable and observed more frequently [48J, so that some further analysis is needed. One 
obstacle has been that practically no examples in higher dimensions are known, and certainly 
not in an explicit fashion. This, in turn, is where the classic Thue-Morse system excels: Since 
its original (spectral) discussion in [47} 136] and its later reformulation in [28J, it has been the 
paradigm of singular continuous (dynamical) spectrum, and it is actually rather easy to also 
calculate its distribution function explicitly and with high precision; see [8] and references 
therein for details and [19] for an approach to the analysis of its 'fractal' aspects. The good 
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accessibility of the latter is a consequence of the underlying Riesz product structure; see [SU] 
for background material. 

However, the situation is not as bad as it appears. A substantial step forward was achieved 
in \32\ [33| [M] for lattice substitution systems by identifying modular coincidence as one viable 
and powerful generalisation of Dekking's criterion |19j for pure pointedness in one-dimensional 
substitutions of constant length. The interesting cases are the systems that fail to possess a 
modular coincidence, such as the bijective substitutions of constant length and their higher- 
dimensional generalisations studied in |20l I21| . However, to the best of our knowledge, no 
serious attempt has been made to investigate a genuinely higher-dimensional example with 
singular continuous spectrum explicitly (by which we mean an example that cannot be written 
as a product of one-dimensional systems, such as those appearing in |26j). 

Here, we report on a planar example that was identified as a good candidate in [22J, and 
could be tackled with the methods explained in [23]. It it known as the squiral tiling and 
appears in [25\ Fig. 10.1.4], where it was constructed as a simple example for a tiling of the 
plane by one prototile with infinitely many edges (and its mirror image). It is obtained by 
means of the simple inflation rule 




together with the corresponding rule for the opposite chirality of the squiral prototile. It has 
the integer inflation multiplier 3 and defines an aperiodic tiling of the Euclidean plane. The 
term 'aperiodicity' is used in its strong version here, hence meaning that the inflation rule 
defines a unique hull (via the closure of the Z 2 -orbit of a fixed point in the local topology) 
with the property that no element of this hull admits a non-trivial period. A larger patch of 
the tiling is shown in Figure [TJ 

As an application of [23J, Frettloh and Sing checked that this tiling has no modular coin- 
cidence, so cannot be a model set, wherefore this is a candidate for a mixed spectrum. In 
what follows, we reformulate this tiling via a topologically conjugate one (by mutual local 
derivability \15\ [3]) and bring it to the simpler setting of block or lattice substitutions. It 
will then be possible to show that its balanced version (with weights 1 and 1 = — 1 of equal 
frequency) has purely singular continuous diffraction spectrum. As a consequence, the dy- 
namical spectrum will be a mixture of a pure point and a singular continuous part, in line 
with the (implicit) conjecture in |20j. 

Our approach (as indicated already) will be constructive, so that we do not only determine 
the spectral type, but also derive the diffraction measure explicitly. In fact, we can identify 
it as a two-dimensional Riesz product. This means that one can formulate a simple recursion 
for a sequence of 2D continuous distribution functions that converge towards the distribution 
function of the squiral measure, and uniformly so (even though the underlying measures 
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Figure 1. A rectangular patch of the squiral tiling with a 'black and white' 
reflection symmetry in the vertical axis. The displayed patch was obtained 
via two inflation steps of (JT]) from the central rectangular seed (marked), the 
latter comprising four tiles of each chirality. 



are absolutely continuous, and can thus only converge to the squiral measure in the vague 
topology). 

In Section this approach is generalised to primitive and bijective block substitutions 
of constant length over a binary alphabet, in arbitrary dimension d. We prove that these 
systems show singular continuous diffraction if the height lattice is trivial. 



2. Squirals and lattice inflations 

The squiral inflation of Eq. (JTJ) is primitive, with inflation matrix M = (45). The latter 
has eigenvalues 9 = 3 2 and 1, with Perron-Frobenius eigenvector (1, 1), which is in line with 
both chiralities occupying the same area and occurring with equal frequency (via reading it 
as left or as right eigenvector). Moreover, if we add one pseudo- vertex at the middle of the 
long edge of the squiral, the inflation tiling is face to face as well. As such, it defines a strictly 
ergodic dynamical system under the shift action of Z 2 . Note that the special (spiralling) 
vertex always falls on the centre of a square that is formed by 4 squirals of the same chirality. 
This also explains the name as a mixture of 'square' and 'spiral'. It is natural to view the 
squiral tilings as a (local) decoration of the square lattice, where the colour simply codes the 
chirality of the square decoration. This correspondence gives rise to a mutually local rule in 
the sense of pj)l [3] as follows. 
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The first step consists of identifying an induced inflation rule on the coloured square lattice, 
which is simply given by 



(2) 



This defines a block substitution that is bijective and of constant length in the terminology 
of [20]. As a coloured tiling, it is thus defined by a primitive inflation rule (with the same 
inflation matrix M as before) . One quickly checks that the inflation action and the derivation 
rule (now applicable in either direction) commute, so that the block substitution defines a 
hull that is mutually locally derivable (MLD) to the squiral hull. It is clear that the former 
is better suited for our further analysis. 

If one iterates the block substitution starting from a single square with its centre as reference 
point, one obtains a sequence of growing (finite) tilings that converges to a tiling of the plane 
which is a fixed point; see Figure [5] for an illustration. This particular tiling has perfect D±- 
symmetry and defines (via its hull, the closure of its Z 2 -orbit in the local topology) a minimal 
dynamical system that is uniquely ergodic. This follows from standard arguments because 
the inflation is primitive and leads to a face to face tiling of finite local complexity. As a 
consequence, the fixed point is linearly repetitive, which implies the claim by |31|, Thm. 6.1]. 
For some general results on tiling dynamical systems, we refer to |41 | I44| , [20]. 

3. Autocorrelation coefficients 

Let us now represent each (coloured) square by a point at its lower left corner that carries 
a matching colour, which we call 1 (for white) or 1 (for grey). The hull X is then a closed 
subspace of {1, 1}^ 2 that is invariant under the shift action of Z 2 by construction. Given 
w G X, we define the corresponding Dirac comb 

(3) u = w5 z2 = w z$z, 

where we write z G 1? as a pair of integers from now on. For the weights, we use the 
convention 1 = — 1. Following the approach pioneered by Hof |27| . we define the natural 
autocorrelation measure 7 of u; as 



(4) 



7 = uj © U) 



wL * wL 
li m ICV 

N^oo (2N + 1) 2 



where Cjv denotes the closed centred square of side length 2N. Here, ]x denotes the measure 
defined by j2(g) = n(g) for g € C C (R 2 ), with g(x) := g{x). The limit in Eq. (j3J) always exists 
due to unique ergodicity, and is independent of the choice of w € X. The measure 7 is of the 
form 7 = T]5 7 2 with coefficients 



7j{m,n) 



lim - -7 

JV-xx) (2N + 1)' 



N 

E 

kd=-N 



w kl w k-rnl-m 
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Figure 2. A square-shaped patch of the squiral tiling with full /^-symmetry, 
as obtained from the equivalent block substitution ([2]). 



where the result is again independent of the choice of w. We discuss the convergence of these 
sums in more detail later. 

Due to the symmetry of our system, it suffices to consider the positive quadrant and 
formulate the autocorrelation accordingly. This simplifies our further calculations. To this 
end, we consider the square inflation with the lower left corner as new reference point, and 
observe that this, when starting from a single block, leads an iteration sequence which fills 
the positive quadrant only. To obtain a complete tiling, one may start from the legal seed 
J i with reference point in its centre. Recall that a patch is called legal when it occurs in 
the n-fold substitution of a single letter, for some n G N. The iteration of the substitution 
then converges towards a 2-cycle that covers Z 2 , each element of which defines the same hull. 
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The substitution is still a 2-cycle when restricted to the positive quadrant. Due to the D±- 
symmetry of the pattern in Figure [21 and hence that of the entire hull X, it is sufficient to 
calculate the autocorrelation coefficients within the positive quadrant. 

Denote the two configurations in the positive quadrant (written in the alphabet {1,1}, 
with 1 = 1) by v and gv. They satisfy 



(5) (gv) 



v m ^ n , if r = s = mod 2, 



\v mn , otherwise, 

where m, n > and < r, s < 2. The autocorrelation coefficients clearly satisfy 

, N-l ^ N-l 

(6) 77(777, n) = Yim^ ^ (e v )k+m,t+n = V M v k+m,e+n 

00 k,e=o 00 k,e=o 

for m,n > 0. All limits exist due to unique ergodicity and the fact that the sum is an orbit 
average of a continuous function Moreover, we have used the symmetry to write rj via 
the positive quadrant only. Clearly, we have 77(0,0) = 1 together with 

(7) rj(—m,n) = 77(777,, —n) = rj(—m,—n) = 77(771,77,), 

which specifies the function 77 on all of Z 2 . Note that the L^-symmetry of our system also 
implies the relation 77(777,77) = 77(77,777). 

For our further analysis, we introduce the shorthand 

(777,77) := 77(777,77). 

Considering the coefficient (3?77+r, 3n+s) for gv and fixed 777, 77 > and < r, s < 2, one can 
split the defining sum modulo 3 and use Eqs. ([5]) and ([6]) to derive the recursion relations 

(3m, 3n) = (m, n) 

(3m,377 + l) = — 1(?77, 77) + |(m, 77 + 1) 

(3m,3n+2) = |(m, 77) — |(m, n+1) 
(3m + l,3n) = — |(tt7, 77) + |(m + l, 77) 

(8) (3777 + 1,377 + 1) = -|((777+ 1, 77) + (777 , 77+ 1)) + i(777 + l, 77 + 1) 
(3777 + 1,377 + 2) = -|((777 , 77) + (777 + 1,77+1)) + i(777 + l, 77,) 

(3777 + 2,377) = I (777 , 77) — I (777+I, 77) 
(3777 + 2, 377 + 1) = -|((777 , 77) + (777 + 1,77+ 1)) + i(777 , 77 + 1) 
(3777 + 2,377 + 2) = i(777 , 77) - |((777 + 1, 77,) + (777 , 77 + 1)). 

All relations are linear and of the form 

min(l,r) min(l,s) 

(9) (3?77+r,377+s) = a^f{m+k,n+£), 

k=o e=o 

where the coefficients are elements of{— 1,0, i, i, 1}. 
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Lemma 1. The autocorrelation coefficients (m,n) of the block substitution system exist for 
all m,n € Z and are completely determined by (0, 0} = 1 together with the linear recursions 
(JSJ), which hold for all m,n € Z. 

In particular, one has the special values (0, ±1) = (±1,0) = — |, (±1,±1) = g, (0, ±2) = 
(±2,0) = (±2, ±2) = ^ and (±1, ±2) = (±2,±1) = Moreover, the coefficients 

satisfy ( — l) m+n {m,n) > for all m,n € Z. 

Proof The existence of the coefficients is clear by unique ergodicity, as discussed earlier. The 
recursion relations (|8|), for m,n > 0, follow from an elementary (though somewhat tedious) 
calculation as mentioned above. The initial condition (0, 0) = 1 is clear, while the other special 
values for non-negative arguments can then be successively calculated from the recursions by 
solving linear equations. As is easily seen, the coefficients (m, n) with m, n > 2 are then 
determined recursively. 

One can now explicitly check that the same method (formally) also gives the other special 
values, which are in agreement with the symmetry relations ((7J) and (m,n) = (n,m). Indeed, 
one can verify that the recursions (extended to all m, n € Z) respect all ^-symmetries. The 
linear recursion relations thus determine (m, n) from (0,0) on the entire lattice Z 2 . 

The positivity claim is certainly true for all —2 < m,n < 2. Inspecting the recursion 
relations (jHj), one sees that all terms on the right hand side of any single relation have the 
same sign, and that the parity on the left hand side matches the formula, so that our claim 
follows inductively. In particular, due to the recursive structure, no coefficient can vanish. □ 

Let us formulate another, rather surprising consequence of the recursive structure, which 
will significantly simplify one of our later estimates. 

Lemma 2. One has (m, 0) 2 — (m,n) 2 > for all m,n S Z. 

Proof. Our claim follows if we show that \{m, n)\ < |(m, 0)| holds for all m,n > 0. Since 
\(m,n) \ = (—l) m+n (m,n) by Lemma [fl the absolute values satisfy the recursion relations ([8]) 
with each coefficient on the right-hand sides replaced by its modulus. The inequalities are true 
for < m,n < 2 (by inspection of the initial condition and the special values of Lemma [1]). 
The claim now follows by induction, where one applies the recursions once to differences of 
the form \(3m+r, 0)| — |(3m+r, 3n+s)\, with < r, s < 2, which are all non-negative. □ 

The autocorrelation coefficients define a positive definite function on Z 2 . By the Herglotz- 
Bochner theorem, compare LemmaOin the appendix, it is thus the (inverse) Fourier transform 
of a unique positive measure on the 2-torus T 2 = [0, l) 2 , so that 

(10) (m,n) = I / exp(27ri(ma; + ny)) dfi(x, y) = / cos(2-7r(mx + ny)) dfj,(x, y). 
Jo Jo Jo Jo 

The second equality follows from the L>4-symmetry of the coefficients, which implies the 

corresponding symmetry for fi. As (0, 0) = 1, we see that /i is a probability measure on T 2 . 

Remark 1. The Herglotz-Bochner theorem can be used to construct examples of positive 
definite functions on Z 2 with Z^-symmetry that do not satisfy the inequalities of Lemma [2j 
For instance, a + (1 — a)(— l) m+n is the Fourier transform of a^ o) + — a )^i i)i but 
violates the inequalities for a G (0, 1) and m + n odd. 
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Before we continue our analysis of the planar case, let us look at an important one- 
dimensional subsystem. It will reappear later in an essential way. Moreover, it serves to 
introduce the methods we use. 



4. A RANK 1 SUBSYSTEM 

Consider e(m) = (m, 0} for m £ Z, which defines a positive definite function on Z, with 
(— l) m e(m) > for all m £ Z by Lemma [2 One has e(0) = 1 together with the recursions 



e(3m) = e(m) 

2 , 1 
— i 
3 

1 , , 2 



'11) e(3m + l) = --e(m) + -e(m + 1) 



e(3m + 2) = -e(m) e(m + l), 

3 9 

which hold for all m 6 Z. We consider the measure 7 = e<5 z on the lattice Z, which corresponds 
to the autocorrelation of our planar system along the x direction (or, by symmetry, along the 
y direction). Its Fourier transform, by [4 S Thm. 1], is of the form 7 = v * 5% with v = 7|[o,i)- 
Following the approach of [28} [6], one can easily see that v is a purely singular continuous 
measure, and it is possible to compute the corresponding distribution function explicitly. 

The absence of a point part of v follows by Wiener's lemma. We provide a proof of the latter 
in the Appendix that is tailored to our later needs; compare [291 Cor. 7.11] or [371 Sec. 4.6]. 
Define Ei(iV) = ^~Je(m) 2 , and recall that |e(n)| < e(0) from positive definiteness. Now, 
using Jensen's inequality via (a + b) 2 < 2(a 2 + b 2 ), we can estimate 

2 N-l 



Ei(3iV) = ^^e(3m + r) J 



r=0 m=0 

N—1 

(12) = Sl (iV)+^(-^(m) + i6(m+l)) 2 +(i e (m)-^(m + l)) 2 

m=0 ° ° 

< (1 + |) Sx(iV) + f^N) 2 - e(0) 2 ) < ^ ^(N) , 

which shows that Ei(iV) = 0(A^ a ) with a = log 3 (133/81) < 1/2, so the measure z/ has 
no points. The absence of an absolutely continuous component can be shown in complete 
analogy to the Thue-Morse case by invoking the Riemann-Lebesgue lemma, so we conclude 
that v is purely singular continuous. The corresponding distribution function «P is defined via 
<P(x) = u([0, x]j for x € [0, 1] and extended to R by means of <&{x + n) = n + <&(x) for n S Z, 
which gives a continuous, non-decreasing function on R. The general properties of <P follow 
by the arguments used for the (generalised) Thue-Morse system; compare [9j [6] for details. 
Let us summarise some of them as follows. 

The continuous distribution function ^ possesses the Fourier series representation 



(13) <P(x) = x+ ^2 ~ — ~ sin(2-7rm2;) = x e(m) sinc(27rmx) , 



rmr 

m>l 
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with sinc(z) = sm ^ and sinc(O) = 1. The first sum is uniformly (but not absolutely) con- 
vergent. The second equality is a consequence of <£> being the integral of a (formal) Fourier- 
Stieltjes series with coefficients e(m); compare [39J, Sec. 1.2.6] and [291 Sec. 7]. The limit can 
be approximated by distribution functions <£ n of absolutely continuous positive measures, 

3«-l , s. 

EC ( TT~l/ ) 
nV ' sm(2irmx) , 
mir 

m=l 

where <&q{x) = x (which is the distribution function of Lebesgue measure A) and the coeffi- 
cients satisfy the initial conditions c (m) = 5 m together with the recursions 

c n (3m) = c n _ x [m) 
c n (3m+l) = -gC n _ 1 (m) + -c n _ 1 (m+l) 

c n (3m+2) = -c^im) - ^^(m+l) . 

The distribution functions <P„ define a sequence of absolutely continuous measures that con- 
verge vaguely to our singular continuous measure defined by Here, the individual Fourier 
series for ^ n (x) — x are finite sums, hence trigonometric polynomials with fundamental period 
1. Moreover, we have uniform convergence <P n — > as n — > oo by the same argument as in 
[BJ, which is based on the stepping-stone argument from [16, Thm. 30.13]. 

The corresponding densities (f n (x), defined by d<P n (x) = <p n (x) dx, have the Riesz product 
representation 

(is) ip n [x) = \{e{^x), 

1=0 

where 6{x) = 1 — | cos(27rx) + | cos(47T2;) is a strictly positive, 1-periodic function that satisfies 
9(1 — x) = 6(x) and 6(x)dx = 1. The functions <P A and (p A are illustrated in Figure EJ 
The measure defined by <P is represented as the infinite Riesz product YY^LqO^x), to be 
understood with convergence in the vague topology. In summary: 

Proposition 1. The function e: Z — > K defined by Eq. (jlip together with e(0) = 1 is positive 
definite. It is the Fourier transform of a probability measure v on [0, 1) that is purely singular 
continuous, with distribution function <P of Eq. f 1 1 3 [) and a representation as an infinite Riesz 
product. □ 

5. Analysis of the planar case 

Since all ingredients to the above arguments are also available in higher dimensions, we 
proceed constructively. Here, with the coefficients from Lemma (H we define 

N-l 

(16) E(iV) := ^ (m,n) 2 , 

m,n=0 

which is now a sum over iV 2 terms. 

Lemma 3. The sums of Eq. {E} satisfy Z(3N) < ^ S(JV) < 4E(iV). 
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Figure 3. Distribution function <&a{x) (left) and its smooth density (pi(x) 
(right), according to Eqs. (fill) and (fl~5l) . 



Proof. The recursion relations ([S]) imply that 
3N-X N-X 2 

m,n=0 m,n=0 r,s=0 

1 

= — Y 84(m,n) 2 + 14(m,n+l) 2 + 14(m+l,n) 2 + 7(m+l,n+l) 2 

m,n=0 

+ 2(4(m,n) - 2(m,n + l) - 2(m+l,n) + (m + 1, n+1}) 2 . 

One can now use Jensen's inequality in the form (a + b + c + (f) 2 < 4(a 2 + 6 2 + c 2 + (i 2 ) and 
observe the relation 

N-X N-X 

Y (m,n+l) 2 = S(iV) + ^((m,7V) 2 - (m,0) 2 ) < E(iV), 

m,n=0 m=0 
which follows from LemmaO The analogous estimate holds for ~L (m + l, n} 2 , while 



AT-l 



(m+l,n+l) 2 = E(iV~) + {(N,N) 2 - (0,0) 2 ) 



m,n=0 



JV-1 



N-X 



+ Y({m,N) 2 - (m,0) 2 ) + Y((N,n) 2 - (0,nf) < E(N). 

m=X n=X 

Putting everything together leads to the estimate 



^1 Q 

E(3AT) < ^ E(iV) 



which implies our claim. 



□ 
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Remark 2. Let us note that, without Lemma [2j one would obtain 

31 Q 

(17) £(3A0 < _S(iV) + 0(iV), 

which is weaker but still sufficient for the application of Wiener's lemma later on. This is the 
type of relation that one can expect in similar examples that fail to satisfy Lemma [2j 

Lemma 4. Let £ : Z 2 — > R be a function that satisfies the recursion relations §8§, with 
C(0)0) — 0. Then, Q is a positive definite function on Z 2 , and defines a unique positive 
measure fx^ on the 2-torus T 2 , via 

C(m,n) = [ e^d^x), 

where /i^(T 2 ) = £(0,0). Moreover, the measure fi^ is absolutely continuous, relative to the 
Haar measure on T 2 , if and only if £(0, 0) = 0. In the latter case, Q = 0. 

Proof. When C(0>0) = 1; we have ((m,n) = (m,n) by Lemma [U which is positive definite 
by construction. Since the recursion © is linear, £(0, 0) = a > leads to £( m i n) = a {m, n), 
which is still positive definite on Z 2 . The Herglotz-Bochner theorem |42j results in the rep- 
resentation via the unique positive measure /i^; compare Lemma [7] below. It is a probability 
measure if and only if £(0, 0) = 1. 

Whenever a > 0, the special values Q{m, n) with —2 < m,n < 2 are different from 0, again 
by Lemma [H wherefore £(3m,3n) = Q{m,n) transports them all the way to infinity. By 
the Riemann-Lebesgue lemma [39} |4"2] , this implies that (j.^ cannot be absolutely continuous 
relative to Lebesgue measure (which is the Haar measure on T 2 ). The only exception is 
C(0, 0) = 0, which forces also all special values to vanish, and hence £ itself. □ 

Theorem 1. The diffraction measure 7 = %j of the balanced Dirac comb oj of Eq. ^ is a 
translation bounded, positive measure that is purely singular continuous. All elements of the 
hull of oj possess the same autocorrelation and the same diffraction measure. 

Proof. The diffraction measure is 7 = fj,*5^ 2 , by an application of [4, Thm. 1], where \i is the 
positive measure from Eq. (110p . From Lemma[3l we know that -^-X(iV) — > as N — > 00, so 
that Wiener's lemma (see Lemma [7] in the appendix) implies that [i is continuous. 

We may now employ the unique decomposition /i = /i sc + fj, ac as a sum of non-negative 
measures, relative to Lebesgue measure on T 2 . Defining the functions r/ ac and rj sc as the 
inverse Fourier transforms of /i ac and /u sc , in analogy to Eq. (|10p . one finds 

r] ac (m, n) + 7] sc (m, n) = (m, n) 

for all (m, n) £ Z 2 , together with ?? ac (0,0) > and ?7 SC (0, 0) > 0. Since fi ac and fi sc are 
mutually orthogonal in the measure sense, it is clear that both rj ac and r/ sc satisfy the same 
set of recursions, namely those of Eq. (|8|), but possibly with different initial conditions. By 
an application of Lemma [U we see that r] ac must vanish, so that (m,n) = r/ sc (m,n) on Z 2 . 
This implies \x, and hence 7, to be purely singular continuous. 

The hull of u), which is its orbit closure in the vague topology, is strictly ergodic. Conse- 
quently, the autocorrelation coefficients of each Dirac comb in this hull are the ones specified 
by Eq. ([8]) and Lemma [U which implies the last claim. □ 
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To further investigate the probability measure fi, and then also the diffraction measure 
7 = n * 5j2i let us define a multi-dimensional distribution function F, compare [24, §23], as 
F(x, y) := /x([0, x) x [0, yu , which is continuous from the left and from below by construction; 
see [TBI Sec. 1.6] or [T71 Sec. 12] for an alternative, equivalent approach. Note that F is non- 
decreasing along any line parallel to the x-axis or the y-axis. The (continuous) measure v 
from Proposition Q] is the marginal of fx in the sense that 

(18) v([0,x)) = n([0,x) x [0,1)) and M ([0,l) x [0,y)) = u([0,y)) 

holds for all x,y £ [0, 1), where the second relation follows by the symmetry of our system. 
Since the distribution function of v is continuous by Proposition [IJ an application of |45} 
Thm. 2.3] shows that F is a continuous function. This is one of the somewhat subtle points 
to observe when using higher-dimensional distribution functions. 

To continue, it is advantageous to extend F to a continuous function on R 2 , which is most 
easily done via F(x,y) = 7([0, x) x [0, y)) for x,y > 0. Here, F is clearly continuous, by an 
extension of the above argument across the lines {x = m} and {y = n}, where F(m, y) = ^{y) 
and F(x,n) = <&(x). In particular, we thus have 

F(x,y) = 7([0,x] x [0,y]) 

for x, y > 0. This can now consistently be extended to all of R 2 by setting F(—x,y) = 
F(x,—y) = —F(x,y) and hence F(—x,—y) = F(x,y). In particular, one has F(0,0) = as 
well as F(0,y) = F(x,0) = 0, and F is continuous on R 2 . 

Lemma 5. Let F be the continuous distribution function defined by F(x, y) = 7([0, x] x [0, yf) 
for x,y > together with F(—x,y) = F(x, —y) = —F(x,y). Then, F satisfies 

F(x + 1, y) = F(x, y) + ${y) and F(x, y + 1) = F(x, y) + $(x) 

for arbitrary x, y € R, where $ is the function from Eq. (|13p and Proposition^ The function 
h defined by h(x,y) = F(x,y) — x$(y) — y<P(x) + xy is continuous and Z 2 -periodic. 

Proof. Observe first that Eq. (|18p has a natural extension the periodic measures v * 5% and 
7 = M*£ Z 2- 

Now, for x, y > 0, one has 

F(x,y + 1)-F(x,y) = / = ' tfy = = d{u * S z ) = 

io io JO JO JO 

where {x} denotes the fractional part of x, and the 1-periodicity of 7 (in y) was used twice, 
in the second step via [{y}, {y} + 1] = [{y}, 1] U (1, {y} + 1] together with the fact that 7 is a 
continuous measure. The other identity follows analogously, and the extension to all of R 2 is 
clear by symmetry. 

The final claim now follows from a simple calculation. □ 

Due to the underlying symmetry, h has a Fourier series of the form 

00 

h(x,y) = a m,n sin(27rmx) sin(27rny), 

m,n=l 
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where a m n = ^i^ n follows from a routine calculation; see [I] for background on multiple 
Fourier series. Together with Eq. fjl3[) . this leads to the series representation 

oo , > 

— sin(27rm:r) sin(27rny) 

m,n=l 

(0> n ) (m,0) 

(19) = xy + x > sm(2-7rny) + y > sm(27rmx) 

Trn ' Tim 

n=l m=l 

oo , , 

E(m, n) 
— ^ sm(27rmx) sm(27rny) , 
7r z mn 

m,n=l 

which is actually uniformly convergent (for summation over square-shaped regions). In line 
with Eq. (|13p . one can also rewrite the distribution function as 

F(x,y) = xy (m, n) sinc(27rma;) sinc(27rny), 



which highlights the structure as an integral over a planar Fourier-Stieltjes series. 

It is important to observe that F is well-defined once the measure \x is given, and that F 
specifies the measure uniquely; compare [T71 Thm. 12.5] for details. We can thus employ the 
Lebesgue-Stieltjes approach to measures also in this more general situation. 

For actual calculations, it is advantageous to use an approximation to F via a sequence 
of distribution functions with densities, pretty much as in the Thue-Morse example. Indeed, 
the recursion relations (JSj) , via an explicit but somewhat tedious calculation, implies the 
functional relation 

(20) F(x,y) = -j J tf(|,|)dF(*,y) 

(written in Lebesgue-Stieltjes notation) with the kernel function 

1 9 

(21) $(x,y) = -(l + 2cos(2vrx) + 2cos(2vry) -2cos(27r(a;+y)) -2cos(27r(x-y))) , 

which is shown in FigureUl Clearly, t? is a non-negative, Z 2 -periodic function that is symmetric 
in both arguments. Moreover, it satisfies the relations — x,y) = y) = 1— y) = 
$(y, x) together with the normalisation Jj ^ 2 $(x, y) dx dy = 1. A general formula for i? will 
be discussed below in Eq. ([30]) . 

The functional relation (|20|) can now be employed to define an iterative calculation of F as 
follows. One starts from F^°\x,y) = xy (which corresponds to the measure dF(°) = A) and 
continues with the iteration 

(22) F^\x,y) = 1 -J^ J tf(|,|)dFW(*, y ). 



14 



MICHAEL BAAKE AND UWE GRIMM 




Figure 4. The function of Eq. on [0, l] 2 . 
The functions have the form 

3 iv_l (JV) 3 iv_i (AT) 

F^ N ^(x,y) = xy + x — — sin(27rny) + y \^ m '° sin(27rrw:r) 

n=l m=l 
3^1 ~(A0 



+ m,n sin(27rmx) sin(27rnj/), 



Trunin 

m,n=l 



where the coefficients f3m,n are defined for N, m, n > by the initial conditions /?m,n = S m0 5 n0 



together with the recursion 

B {N+i) V- V « (r,s) /5 



min(l,r) min(l,s) 

(iv+i) _ fl W 

fc=0 £=0 



with JV > and the same coefficients oijif' as in Eq. ([9]). This can once again be verified by 
a direct calculation. 

All F^ N > represent absolutely continuous measures, so that we can define Radon-Nikodym 
densities via dF^ N \x,y) = f( N \x, y) dxdy, where 

ax ay 

This gives = with the function i? from Eq. (|2ip . via the application of some trigono- 
metric identities. The iteration now reveals the Riesz product formula 

N-l 

(23) f (N) (x,y) = HH^x,^), 
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Figure 5. The distribution function of Eq. ([22]) (left) and the corre- 
sponding density of Eq. (|23|) (right), approximating the diffraction mea- 
sure of the squiral tiling (in the version of Figure [2]) on [0, l] 2 . 

which highlights the deeper role of t? in Eq. (|22p . The distribution function and the 
corresponding density are shown in Figure [5j 

Remark 3. The measure defined by the continuous function F is represented by the infinite 
Riesz product Yle^o x > 3 f y), to be understood with convergence in the vague topology. By 
a comparison of Eqs. (j20[) and (122|) . it is clear that F is a fixed point of the latter. In fact, 
within the class of distribution functions with certain continuity and additivity constraints, 
it is the only fixed point, with uniform (but not absolute) convergence towards it from the 
initial condition F(°\ again by arguments analogous to [16, Thm. 30.13]. 

6. A TOPOLOGICAL FACTOR WITH MAXIMAL PURE POINT SPECTRUM 

It is well known that the one-dimensional Thue-Morse system admits a topological factor 
with maximal pure point spectrum, which can be defined by the period doubling substitution. 
In fact, the latter is induced by a simple sliding block map of width 2. The corresponding fac- 
tor map is globally 2-to-l between the hulls; see [6] for details and an extension to generalised 
Thue-Morse sequences. Amazingly, a similar approach also works for the squiral tiling. 

Define the mapping ip: {1, l} 1,2 — > {1, l} z2 by w i— > %{jw with 

which is continuous. It is clear that Y := defines a factor for the action of 1? . Moreover, 
the inflation rule (pQ) induces a new inflation on Y, which reads 

1 1 a 

(24) 111 

a i — > 1 1 1 

with a € {1,1}- Note that the arrangement matches that of the inflation rule ([I]). One 
checks consistency by starting from the 14 legal patches with 2x2 squares, and verifies that 
they inflate to larger patches that produce the correct 3x3 blocks in the lower left corner 
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under ip. One can check that the block map is not globally 2-to-l, though it is 2-to-l almost 
everywhere. For a more detailed discussion and a classification of the substitution factors of 
the squiral, we refer to [7]. 

Consider one of the two possible fixed points of the induced substitution (|24p and decompose 
1? = A + OA^, where A± = {(m,n) £ 1? \ w m ^ n = ±1}. Using e 1 = (1,0) and e 2 = (0,1), 
the fixed point property now induces the set- valued relations 

(25) A± = (3Z 2 + 5±) U (3A± + 2(e x +e 2 )) 

with the translation sets S+ = {2e 1 , 2e 2 , 2e 1 +e 2 , e 1 +2e 2 } and S- = {0, e x , e 2 , e 1 +e 2 }. 
The solutions are 

4b = (J ( 3 " z2 + 3 "" l5 ± + ^ ~ 1 )( e i+ e 2)) u A± 

n>l 

where either = {— (e 1 + e 2 )} and A_ = 0, or vice versa. The special role of the point 
(— 1, —1) becomes transparent in the topology of the 3-adic numbers. 

The sets A± are model sets with two copies of the 3-adic numbers as internal space. The 
Dirac combs 5 A± are thus pure point diffractive by the model set theorem; see [4~3| Thm. 4.5] 
or [14^ Thm. 2]. A standard calculation now reveals that the corresponding Fourier module 
is Z[^] x Z[|]. By general arguments [SHU], this equals the dynamical spectrum of Y. 

By the results of |20| . the pure point part of the dynamical spectrum of X thus equals the 
dynamical spectrum of Y, which shows that Y is indeed a topological factor that exhausts 
the pure point spectrum of the squiral tiling. To summarise, we have shown the following. 

Theorem 2. The squiral tiling dynamical system has a dynamical spectrum of mixed type, 
with pure point part Z[|] x Z[|] and an additional singular continuous part. The latter is 
exploited by the diffraction measure of the balanced Dirac comb of Eq. ([3]) , while the former 
can be recovered from the pure point diffraction spectrum of the factor system defined via 
Eq. (PU>- □ 

7. Generalisations 

When inspecting the structure of the above constructive approach, it is apparent that the 
method is not restricted to the particular example of the squiral tiling or its equivalent block 
substitution rule. Instead, one might as well consider a d-dimensional extension to bijective 
block substitutions on a binary alphabet as follows. 

Let K be a matrix array of dimension K\ x K2 x . . . x Kj, with all Ki > 2 and entries ±1 
(with 1 = — 1 as before), not all equal. Now consider the rule 

q: 1 1 — ► K, 1 1 — >K, 

where K is obtained from K by flipping all signs. This defines a primitive block substitution 
rule with substitution matrix M = (^^), where m is the cardinality of Is in K and fh that 
of Is. The PF eigenvalue is n«^' with (left and right) eigenvector given by (1,1). This is 
in line with a tiling interpretation, with unit cubes of two colours as prototiles. They occur 
with equal frequency in any fixed point tiling. 
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Select any legal block patch of size 2 x 2 x . . . x 2, which contains 2 d blocks, and take its 
centre as reference point (meaning that each sector of % d contains one block of it). This way, 
the origin is a corner of a block. Starting from such a seed, we can now iterate g to create 
a sequence of block patches that ultimately cover all positions of Z d , due to our assumption 
that all Ki > 2. Note that the block patches of type K and K fit together in a face-to- face 
manner in this process. 

Without loss of generality, we may now assume that we have constructed a fixed point 
in this way, which is possibly only true after replacing g by a suitable power of it, and 
modifying K accordingly. This modification is always possible, as follows from an application 
of Dirichlet's pigeon hole principle to the finite set of legal seeds. Note that any power of g 
still satisfies our basic assumptions, wherefore this step is immaterial for the arguments to 
follow (although it might technically complicate explicit calculations due to the increased size 
of the new array K). 

Let w be the fixed point under consideration. It is then specified by the recursion relations 

( 26 ) w K 1 m 1 +r 1 ,...,K d m d +r d = K r 1 ,...,r d w m 1 ,...,m d 

where all rrn € Z and < r% < Ki, while the coefficients n r r € {±1} are the elements of 

1 '***' d 

K, numbered accordingly. Due to the primitivity of g, the fixed point w defines a repetitive 
configuration in the shift space {±1} Z , wherefore it leads to a hull 

X = X(w) := {t + w | t i Z d } 

that is a minimal subshift of {±1} Z . Moreover, (X, Z d ) defines a strictly ergodic dynamical 
system by standard arguments [4lJ HU [20J HO] . 

If u> = w5 zd is the attached Dirac comb, its autocorrelation 7 = oj®lo (defined in complete 
analogy to Eq. (|4|)) exists and is of the form 7 = rj5 zd with coefficients 

^ N N 

(27) ri( mi ,...,m d ) = lim^ ——— £ ••• £ ^^w^^^. 

V ' h=-N k d =-N 

The sums are orbit averages of a continuous function under the action of Z d , which always 
exist due to unique ergodicity (by an application of the stronger version of Birkhoff 's ergodic 
theorem |46j). 

One can now repeat the constructive calculation of the squiral block substitution in this 
more general setting. 

Lemma 6. Let g be a primitive, bijective block substitution on a binary alphabet, with a fixed 
point configuration w that satisfies Eq. (I26h . The corresponding autocorrelation coefficients 
satisfy rj(0, ... ,0) = 1 together with the linear recursion relations 

Ti(K 1 m 1 + r 1 ,...,K i m d +r i ) = V • • • ( l! '"'^ d ) r){m x + s x ,. . . ,m d + s d ) 

s 1= s d =0 \ S l'---' S d/ 

with rational coefficients (5^!!!^) of modulus < 1. This recursion determines all coefficients 
uniquely once the initial condition 77 (0, ... , 0) is given. 
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In particular, one has ( Sl ,'.'.'.'s (J ) = o ' • • • ' ^s d o> so that the relation 

■n{K x m x ,...,K d m d ) = rj{m x , . . . , m d ) 
holds for all (m x , . . . , m d ) £ Z d . 

Proof. It is clear that r/(0, ... ,0) = 1. For a general r](K 1 m x + r 1 , . . . ,K d m d + r d ), we start 
from the right-hand side of Eq. (127p (for finite N) and split the summation over the ith index 
modulo Ki. For each term, we can now apply the recursion (|26[) . The sum can then be 
regrouped to contribute to coefficients of the form rj(m x +s x , . . . ,Tn d +s d ) with < s i < 1. In 
the limit as N — > oo, the contributions sum up to 

/ \ ^1-1-^(1-^) K d -l-r d (l~s d ) 

( , f ri ,...,r d \ ^ ^ K t 1 ,...,t d \+r 1 -s 1 K 1 ,...,t d +r d -8 d K d 

which follows by carefully keeping track of the various overflow rules due to the calculations 
modulo K{. The claim on the rationality and boundedness is then obvious. 

Inspecting the equations for the autocorrelation coefficients, one notices that for some of 
them (for instance those with m 1 = . . . = m d = and one T{ = 1) the resulting equation 
is a linear equation for a single coefficient that occurs on both sides. The prefactor of this 
coefficient on the right hand side can never be 1, as one can see from Eq. (|28p together with 
the assumption that the K-array does contain entries of both signs. So, this type of equation 
always has a unique solution in terms of r/(0, . . . , 0), which fixes this particular r\ value. The 
same argument applies to all equations of this type. The remaining coefficients are then 
determined recursively, and the solution is unique as soon as r/(0, ... ,0) is specified. 

The final claim is clear for s x = . . . = s d = 0. If s- = 1 for some i, the corresponding sum 
is empty, and the coefficient vanishes. □ 

By construction, the mapping r/ : Z d — > R is a positive definite function on Z d , with 
?7(0, . . . , 0) = 1. By the Herglotz-Bochner theorem, see [12] or Lemma [7] below, there is then 
a unique probability measure \i on T d = [0, l) d such that 

(29) i](m) = / ••• / exp(27rima:) d/i(x), 

Jo Jo 

where we write m = (m 1 , . . . , m d ) etc. in the remainder of this section, with mx denoting the 
scalar product. The positive measure \i has a unique decomposition into three components, 

V = A^pp + /^sc + ^ac> 

with Lebesgue measure (as the Haar measure on T d ) as reference measure. 

Proposition 2. Let g be as in Lemma® with autocorrelation coefficients r\ and associated 
positive measure \x onT d . Then, the spectral type of [i is pure, which means that it is either a 
pure point measure, a purely singular continuous measure, or a purely absolutely continuous 
measure. 
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Proof. By standard arguments, we know that n pp _L ^ sc _L fi ac ± fi pp in the measure-theoretic 
sense. Now, the recursion relations imply a set of functional relations for ^u, which must hold 
for each component separately as a consequence of the orthogonality relations. 

Defining r] a (m) = fi a (m) for a € {pp,sc,ac} separately, each type of autocorrelation 
coefficient must then satisfy the same recursion relations as r] itself, subject to the condition 

V(0) = >7pp(0) + »7sc(0) + %c(0) = I- 

Note that, due to the linearity of the recursion relations, this splitting indeed leads to a 
solution of the recursion with the correct initial condition, hence the unique solution stated 
in Lemma [6l 

This observation means that the three components of rj{m) are proportional according to 
their individual initial conditions r] a (0), so cannot produce different spectral type. This is 
only compatible with the claim made, and the type is pure. □ 

Let us take a closer look at the spectral types. Assume that the measure \i is absolutely 
continuous. Then, the Riemann-Lebesgue Lemma tells us that rj(m) = r/ ac (m) — > as 
\m\ — > oo. Recalling the last claim of Lemma [6l it follows that we must have rj(m) = for 
all m 7^ 0, and hence r/(m) = 5 m . This would imply \x to be Lebesgue measure on T . A 
careful inspection of the recursion reveals that this outcome is only possible if all coefficients 
of type Co'''''o d ) vanish. However, we have 

(K x -\,...,K A -\\ = K o,...,o K Kl -i,...,K d -i 
V J K l -...-K d T 

which excludes this possibility. Consequently, the measure in Proposition [2] is a singular 
measure. 

Put together, we have proved the following result. 

Theorem 3. Let g be a primitive, bijective block substitution on 7L d , for a binary alphabet. 
Assume that it is extensive in all directions, and let w be a fixed point configuration that 
satisfies Eq. ([26]) onZ d . 

If i] is the corresponding autocorrelation function on % d , and /i the attached positive mea- 
sure on T d according to Eq. flZJJl, the measure \i is singular. In particular, it is either pure 
point or purely singular continuous. □ 

Both possibilities clearly occur. Examples for pure point cases can be constructed via 
substitutions that have fully periodic fixed points. Such fixed points usually have non-trivial 
height lattices, and are thus of limited interest; see \20\ Sec. 3.1] for the precise definition 
of the height lattice of a Z rf lattice substitution. If the height lattice of a primitive bijective 
substitution is trivial, one can construct a cyclic space in the complement of the span of the 
eigenfunctions [20l l2T] . Its spectral measure is nothing but our measure fi, which is then 
inevitably singular continuous. 

Corollary 1. If the block substitution from Theorem [3] has trivial height, the diffraction 
measure of the corresponding balanced Dirac comb is purely singular continuous. □ 
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Remark 4. Singular continuity of a (positive) measure \i comprises rather different situations 
in d > 1 dimensions. For instance, \i can be genuinely singular continuous, such as the measure 
of our squiral tiling above. However, one can also find product measures of a different kind. 
As an example, consider the primitive, bijective block substitution defined by 

1 1 1 
1 i — ►111 

in symbolic notation on the binary alphabet {1,1}. This rule defines a hull that is 2-periodic 
in the horizontal direction and Thue-Morse in the vertical one. Consequently, the measure /U 
in this case is a product measure, written as [i = fi 1 <g> [A 2 , where fi 1 is a pure point measure 
and fj, 2 is the classic TM measure. Nevertheless, fi is purely singular continuous as a measure 
on T 2 . This is an example with a non-trivial height lattice that still shows purely singular 
continuous diffraction, though the corresponding distribution function is not continuous. 

Let us assume that the probability measure \i from Eq. (|29|) is purely singular continuous 
and that the distribution function F for 7 = fi * 5^ d is continuous (which excludes examples 
such as those mentioned in Remark U]). The latter condition is equivalent to the d one- 
dimensional distribution functions F(l, . . . , 1, Xj, 1, . . . , 1) being continuous, by an application 
of \45\ Thm. 2.3]. Then, the recursion relation for the autocorrelation coefficients n can 
meaningfully be turned into an iteration, as in our planar example. One can then derive 
an explicit representation of /i (and hence also 7) as a (generalised) Riesz product in this 
case. In particular, iterating the corresponding functional equation once with initial condition 
F(°\x) = Y\i=i x i leads to F^ l \ with density function 

which plays the role of the trigonometric function $ in this more general setting. It can 
directly be calculated by means of the recursion coefficients as 

(30) £ ••• (n o n cos ( 27r ^)> 

which defines a non-negative function on M. d with Jj d i?(as) dx = 1. This leads to f( N \x) = 
Yle^Q 1 ^(-^i x i) ■ • • i Kd x d) an d thus to an infinite Riesz product for the diffraction measure 7. 
The distribution function can also be written as 

d 

F{x) = x 1 -...-x d ?7(m) sinc(27rmjXj), 
mez d i=l 

in complete analogy to our previous formulas of this type. 

A complete classification of the cases with pure point spectrum and an extension to larger 
alphabets remain as interesting open problems. Another question concerns a better under- 
standing of the topological factors of bijective substitutions. We expect a complete corre- 
spondence between the original dynamical spectrum and the diffraction spectra of the system 
and its factors. 
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Appendix 

Let us formulate Wiener's lemma for the group which is a straight-forward generalisa- 
tion of the case d = 1 in (37J Sec. 4.16]. Here, C n = {u € W 1 \ Wu]]^ < n} denotes the closed 
cube of sidelength 2n, centred at the origin. One has card(Z d n C n ) = (2n+l) d . 

Lemma 7. If rj : Z d — > C is positive definite, there is a unique positive measure on 
T d ~ R d /Z d such that 

rj{k) = [ e 2niku dfi{u) 
Jl d 

holds for all k G Z d . When r/(0) = 1, the measure fi is a probability measure on T . 
Moreover, one has the relation 

v ' kez d nc n teJ d 

where the last sum runs over at most countably many points. 

Proof. Since 7L d and T d form a mutually dual pair of locally compact Abelian groups, the first 
claim is just the Herglotz-Bochner theorem for this situation, compare [32], with fi(T d ) = r](0). 
Now, define g n {u) = ( 2w +i) d Sfcgz d nc n e 2lTlku , which is bounded by 1. It is not difficult 

to see that g n ( u ) n ~ 5 "°°> l{o}(^) holds for any u £ T d = [0, l) d . Using Fubini's theorem and 
noting that p, = /j,, one obtains 



y ' k& d nc n y ' k& d r\C n 

g n {u - v) &n(u) &n(v) w ~ > °°> / l{o}(^~ v) dfi(u) dfi(v) 

T d xT d JT d xT d 

where the last step follows from dominated convergence; compare [38, Thm. 6.1.15] for a 
formulation in sufficient generality. The integral over the diagonal evaluates via Fubini's 
theorem as 

2 



l { o } («-v)d^(«)dAi(t;) = / n({v})dfi(v) = 

T d JT d JT d teJd 

where the last sum (since \i is a positive measure) runs over all points t £ T d with fi({i}) > 0, 
which are (at most) countably many. □ 
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